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ABSTRACT 

It has recently been shown that relaxed spherically symmetric dark matter halos develop from 
the inside out, by permanently adapting their inner structure to the boundary conditions imposed 
by the current accretion rate. Such a growth allows one to infer the typical density profiles of 
halos. Here we follow the same approach to infer the typical spherically averaged profiles of 
the main structural and kinematic properties of triaxial, anisotropic, rotating halos. Specifically, 
we derive their density, spatial velocity dispersion, phase-space density, anisotropy and specific 
angular momentum profiles. The results obtained are in agreement with available data on these 
profiles from iV-body simulations. 

Subject headings: cosmology: theory — dark matter — galaxies: halos 



1. INTRODUCTION 

High-resolution cosmological simulations show 
that relaxed cold dark matter (CDM) halos are 
close to triaxial homologous systems (Bailin & 
Steinmetz 2005) with a variety of axial ratios but 
essentially universal profiles. That is, the shape 
of the spherically averaged radial profile of any 
given structural and kinematic property is always 
very similar, independently of halo mass, epoch, 
environment and cosmology considered; only the 
scaling may depend on such particularities. 



1 Associated with the Institute, de Ciencias del Espacio, 
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Dubinski & Carlberg (1991) and Crone et 
al. (1994) noted that halos of very different masses 
show similar scaled density profiles. Navarro et 
al. (1997, hereafter NFW) then showed that the 
spherically averaged density profile is always well 
fit, down to about one hundredth the virial radius 
R, by the simple expression 

(p)(r)oc 1 (1) 
r(r s + r) z 

with the scale radius r s correlating with the to- 
tal halo mass M within R in such a way that 
the smaller M, the higher the concentration c = 
R/r s , a correlation that was interpreted (NFW; 
Salvador-Sole et al. 1998) as due to the fact that 
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less massive halos typically form earlier when the 
mean cosmic density is higher. 

Although there is nowadays general agreement 
on the previous universal density profile, some au- 
thors claim that higher resolutions yield slightly 
steeper central cusps (Fukushige & Makino 1997, 
2001; Moore et al. 1998, 1999; Ghigna et al. 2000; 
Jing k Suto 2000; Diemand et al. 2004), while oth- 
ers advocate rather the opposite, that the density 
profile becomes shallower as smaller and smaller 
radii are reached (Taylor & Navarro 2001; Power 
et al. 2003; Fukushige et al. 2004; Hayashi et 
al. 2004). Zhao (1996) has proposed a more gen- 
eral practical expression that accounts for all these 
possibilities. Besides, it has more recently been 
shown (Navarro et al. 2004; Merritt et al. 2005; 
Merritt et al. 2006; Graham et al. 2006; see also 
Stoehr et al. 2002) that a new expression of the 
Sersic (1968), in 3D, or Einasto (1965) law 



(P)(r) 



Po exp 



l/r, 



(2) 



with strictly no central cusp, yields still better fits 
to the mass distribution of simulated halos. Down 
to the resolution-limited radii reached by current 
simulations, the discrepancies among all these an- 
alytical fits are of the order of the deviates found 
using any individual one (Dehnen & McLaughlin 
2005; see also Navarro et al. 2004 and Fukushige et 
al. 2004), which explains such a diversity of opin- 
ions. In respect to the c(M) relation, different an- 
alytical expressions have also been proposed that 
try to recover the results of numerical simulations 
at various redshifts (NFW; Eke, Navarro & Stein- 
metz 2001; Bullock et al. 2001a; Zhao et al. 2003). 
Again, there is overall agreement between them 
although the discrepancies become substantial as 
one gets apart from the mass and redshift ranges 
analyzed in simulations. 

On the other hand, the 3D velocity dispersion 
profile a(r) is well fit by the solution of the Jeans 
equation for spherical, isotropic systems resulting 
from the empirical NFW density profile and van- 
ishing velocity dispersion at infinity (Cole & Lacey 
1996; see also Merritt et al. 2006 for the Einasto 
profile). What is more interesting, as noted by 
Taylor & Navarro (2001; see also Ascasibar et 
al. 2004; Rasia et al. 2004; Barnes et al. 2006), 
the (pseudo) phase-space density profile, (p}/a 3 , 



is always close to a pure power law in radius 
(p)(r) 



a 3 (r) 



= Ar~ 



(3) 



with index v w 1.9 (see also Dehnen & McLaughlin 
2005 for a similar relation applying to the radial 
component of the velocity dispersion). 

In addition, Hansen and Moore (2006) have re- 
cently shown that the pressure-supported aniso- 
tropy profile, 



0(r) 
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a 2 (r) 



(4) 



defined as usual in terms of the ID radial and 
tangential velocity dispersion profiles, cr r (r) and 
<7t(r), respectively, is related to the logarithmic 
slope of the spherically averaged density profile 
through the simple linear relation 



dln(p) 
dlnr 



(5) 



with constants a and b respectively equal to « 
-0.2 and w 0.8 (Hansen & Stadel 2006). 

The simplicity of the relations ((3]) and ([5|) might 
suggest that they are more fundamental than the 
universal halo density and velocity dispersion pro- 
files themselves. However, they appear to be 
equivalent: not only do the latter profiles satisfy 
the previous equations, but they are also the only 
profiles to do so. Indeed, the NFW density profile 
is the only physically acceptable, realistic (with 
no central hole), profile with Zhao's (1996) gen- 
eral form that solves the Jeans equation satisfying 
the relation (J3]), both in the simple isotropic case 
(Austin et al. 2005; see also Taylor & Navarro 2001 
and Hansen 2004) and the general anisotropic one, 
provided in this latter case the additional relation 
© (Dehnen & McLaughlin 2005). 

Finally, cosmological simulations show that re- 
laxed CDM halos have a small angular momen- 
tum, typically orientated along the minor axis 
of the inertia ellipsoid, whose modulus J is log- 
normally distributed (e.g., Barnes & Efstathiou 
1987; Ryden 1988; Warren et al. 1992; Cole & 
Lacey 1996; Bullock et al. 2001b, hereafter BOlb; 
Bailin & Steinmetz 2005) in terms of the dimen- 
sionless spin parameters, 



A 



J\E\^ 2 
GM 5 / 2 ' 



(6) 
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or 



V2MRV C Ms/IGMR ' 



(7) 



respectively defined by Peebles (1980) and BOlb, 
where G is the gravitational constant, E is the to- 
tal energy of the halo (with vanishing potential at 
infinity) , R is the virial radius and V c = y/GM/R 
is the circular velocity. According to their re- 
spective definitions, both spin parameters would 
coincide provided halos were singular isothermal 
spheres, but for halos endowed with the NFW den- 
sity profile, one has (BOlb) 



A' w XF(c) , 



with 



0.7 



(8) 



(9) 



The mean A' value appears to be essentially con- 
stant in time (Hetznecker & Burkert 2006) but 
dependent on halo mass, while the mean A value 
is independent of mass (BOlb) but shows a slight 
variation in time (Hetznecker & Burkert 2006) 
due, in principle, to the evolution of halo concen- 
tration and mass distribution. 

The local specific angular momentum vec- 
tor within halos tends to be everywhere aligned 
(BOlb), the cumulative mass-distribution of spe- 
cific angular momenta M(< j) being well fit by a 
simple two-parameter function and the spherically 
averaged specific angular momentum profile (j) (r) 
fairly well fit by the simple expression 



(j)(r) = U)(R) 



M(r) 



M 



(10) 



where M{r) is the mass inside r and index s takes 
values around 1.3. 

As argued by Austin ct al. (2005), the equi- 
librium condition alone cannot explain all these 
universal trends because the Jeans equation is not 
restrictive enough, nor can the initial conditions, 
as very different protohalos lead to essentially the 
same final structural (Austin et al. 2005; Romano- 
Diaz et al. 2006) and kinematic (Hansen & Moore 
2006) properties. Thus, what can only be at their 
origin is the way these systems grow and, in the 
case of the specific angular momentum profile, the 
effects of tidal torques suffered during that process 
(Doroshkevich 1970; White 1984). 



Two extreme points of view have been investi- 
gated. Some authors have looked at the possibil- 
ity that the universal density profile arises essen- 
tially from the effects of repeated mergers (Syer & 
White 1998; Salvador-Sole et al. 1998; Subrama- 
nian ct al. 2000; Dckcl ct al. 2003). Others have 
focused on smooth accretion through spherical in- 
fall (Avila-Reese et al. 1998; Nusser & Sheth 1999; 
del Popolo et al. 2000; Manrique et al. 2003; As- 
casibar et al. 2004). Similarly, the origin of the 
universal specific angular momentum profile has 
been studied in the spherical collapse approxima- 
tion (BOlb) or by considering the cumulative ef- 
fect of mergers (Gardner 2001; Mailer et al. 2002; 
Vitvitska et al. 2002). 

A notable result found along this latter line of 
research is that the density profile predicted from 
spherical accretion resembles greatly that found in 
numerical simulations. However, major mergers 
cannot be ignored in any realistic hierarchical cos- 
mology, which gives little credit to the "pure" ac- 
cretion scenario. Yet, the density profiles of halos 
do not depend on the epoch they suffered the last 
major merger (Wechsler et al. 2002) nor, in gen- 
eral, on their particular aggregation history (Huss 
et al. 1999; Romano-Diaz et al. 2006). Manrique 
et al. (2003) pointed out that this would be well 
understood if the inner structure of relaxed ha- 
los were completely fixed by the boundary con- 
ditions imposed by current accretion and did not 
depend on the halo past history. In fact, taking 
into account that halos develop from the inside out 
during accretion phases (Salvador-Sole et al. 2005; 
Romano-Diaz ct al. 2006; Lu et al. 2006), their 
typical density profile can be readily derived, in 
the spherically symmetric case, from the accretion 
rate characteristic of the cosmology under consid- 
eration (Manrique et al. 2003). 

This accretion-driven density profile is in very 
good agreement with the results of numerical sim- 
ulations in the whole radial, mass and redshift 
ranges reached by them (Salvador-Sole et al. 2007, 
hereafter SMGH) and recovers all the correla- 
tions shown by simulated halos such as the mass- 
concentration relation (Salvador-Sole et al. 2005). 
Furthermore, as recently shown by SMGH, all the 
conditions required in the derivation of this pro- 
file, but the spherical symmetry assumed for sim- 
plicity, emanate directly from the very nature of 
standard CDM. The fact that CDM is collision- 
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less guarantees that the spatial distribution of any 
physical quantity is continuous at any derivative 
order. As a consequence, all halos with mass 
M at the time t undergoing the same accretion 
M(t) during some finite time interval around t 
have identical (non-scaled) M(r) and, hence, p(r) 
profiles regardless of their individual past history. 
On the other hand, the non-decaying, non-self- 
annihilating and dissipationless nature of CDM 
added to the fact that it accretes slowly onto ha- 
los guarantees their inside-out growth. Finally, the 
CDM power-spectrum characteristic of the partic- 
ular cosmology considered fixes the typical accre- 
tion rate and through it the typical density profile 
of relaxed halos of any given mass at any given 
cosmic time. 

Yet, this success would be of limited interest 
if the same explanation did not also hold for the 
mass distribution in more realistic, triaxial halos 
and for any other property apart from the den- 
sity. In the present paper, we show that this 
approach allows one to understand all the uni- 
versal structural and kinematic trends of triax- 
ial, velocity- anisotropic, rotating halos. In Sec- 
tion [21 we show that the total values at a given 
cosmic time of the extensive quantities character- 
izing relaxed halos and the rates at which they in- 
crease during any finite time interval around that 
moment determine completely the corresponding 
spherically averaged profiles. This is used, in Sec- 
tion [3J to infer the typical shape of these profiles. 
Our results are summarized and discussed in Sec- 
tion [H Throughout the present paper we adopt 
the concordance cosmological model characterized 
by {n m ,n A ,h,a s ) = (0.3,0.7,0.7,0.9). 

2. ACCRETION RATES AND UNIQUE- 
NESS OF PROFILES 

Consider a relaxed (i.e. quasi-steady), triaxial, 
velocity-anisotropic, rotating dark matter halo. 
The mass inside radius r takes the usual form for 
spherical systems, 

M(r) = 4tt [ dff 2 {p){f), (11) 
Jo 

in terms of the spherically averaged density profile, 

(p)(r) = hlJ* Lip F* B sin ^M,^), (12) 



where p(r,6,ip) is the local density in spherical co- 
ordinates r, 9 and ip, centered at the peak density 
and with z axis orientated along the total angular 
momentum J within R. 

The kinetic energy inside r also takes the form 
for spherical systems (see the detailed derivation 
in the Appendix) 

K(r) =4tt f r dff 2 (p)(f)^P-, (13) 
Jo 2 

in terms of the usual velocity-non-centered spa- 
tial velocity dispersion profile, <r(r), which coin- 
cides with the spherical average of the local veloc- 
ity dispersion a(r, 6, ip) in the r shelQ. Note that 
in dealing with non-centered velocity dispersions, 
the influence of rotation on the total kinetic en- 
ergy of the system is included only implicitly in 
equation 1X3]) . 

Simulations show that the local angular mo- 
mentum vector is essentially aligned all over the 
halo, eventually except for the innermost and out- 
ermost regions respectively affected by the limited 
resolution and possible boundary effects due to the 
inclusion of infalling matter (BOlb; Bailin & Stein- 
metz 2005). Consequently, the modulus J of the 
total angular momentum vector inside r can be 
written in terms of the modulus j of the local spe- 
cific one as 

J(r) — I dr r 2 dip d6 sin6j(r, 6, <p) p(r, 9, <p) , 
Jo Jo Jo 

(14) 

Thus, using the spherical average profile, 

(i)( r )= , i \ ( \ l °V/ d6 sin8 j(r, 6, cp)p(r, 6, ip), 
^{P)(r)J a Jq 

(15) 

it takes again the same form as for spherically sym- 
metric systems, 

J(r)=A% [ drr 2 (j)(r)(p)(r). (16) 
Jo 

To conclude this preamble, let us add that 
the similarity with the spherically symmetric case 
when dealing with spherically symmetric profiles 



In contrast, the velocity- centered velocity dispersion at r, 
defined as the rms velocity deviation from the here non- 
vanishing mean value, differs in general from the spherical 
average of the corresponding local quantity. 



4 



does not stop here. It also concerns, although only 
approximately, other profiles such as the poten- 
tial energy inside r or relations such as the Jeans 
equation and the scalar virial relation (see the Ap- 
pendix). In particular, from equations (|A10[) and 
(|A31|) one has the relation 



1 



R 3 (p)(R)Jo 



dr r 2 (p)(r) 



GM(r) 



(17) 



2.1. Density Profile 



As mentioned, relaxed halos evolve inside-out 
between major mergers, that is, they keep their in- 
ner density distribution unaltered and just stretch 
it outwards (Salvador-Sole et al. 2005; Romano- 
Diaz et al. 2006; Lu et al. 2006). This is the natu- 
ral consequence of the fact that centered spheres of 
arbitrary radii conserve the mass (standard CDM 
is non-decaying and non-self-annihilating) and en- 
ergy (standard CDM is dissipationless) and that 
the characteristic accretion time of halos is sub- 
stantially smaller than their dynamical time (see 
SMGH). It is true that relaxed halos are triaxial, 
so they can in principle suffer tidal torques from 
the surrounding matter affecting the kinetic en- 
ergy of such inner spheres. However, the tidal field 
produced by the surrounding anisotropic large 
scale mass distribution is very stable (possibly 
except for the short time interval before major 
mergers) and relaxed halos remain elongated along 
one fixed privileged direction. Consequently, tidal 
torques have a minimal effect on the internal kine- 
matics of relaxed halos; they are only important 
near maximum expansion of protohalos where the 
angular momentum of those structures is gener- 
ated. 

The inside-out growth implies that, at any time 
t during such accretion periods, the total mass is 
given by (see eq. ([TT] ) 



M(t) = M[R(t)] = 4tt 



R{t) 



drr 2 (p)(r), (18) 



where the function (p) inside the integral on the 
right is independent of time. Thus, by differenti- 
ating eq. (HHJ) we are led to 



(p)(t) = (p)[R(t)] 



M 



4irR 2 (t)R 



(19) 



As the virial radius encompasses a region with in- 
ner mean density equal to some factor A v i r (t) (e.g. 
Bryan and Norman 1998) times the mean cosmic 
density pit), 



R(t) = 



3M(t) 



"I 1/3 



47rA vir (t)p(t) 



(20) 



R in equation ([19]) is a function of M. Then, 
all relaxed halos with the same value of M (and 
R) at t, accreting mass at the same rate M dur- 
ing any finite time interval around t, will develop 
identical (p) profiles over the corresponding finite 
radial range. As CDM is collisionless and free- 
streaming and, hence, it admits no discontinuity 
in the steady spatial distribution of any structural 
or kinematic property, the function (p) (r) must 
be analytical. Consequently, if the density pro- 
files of such halos coincide over some finite radial 
range, they necessarily do at any other radius. As 
explained in SMGH, this means that the density 
profile of relaxed halos permanently adapts to cur- 
rent accretion and, hence, does not depend on the 
halo's past aggregation history. 

2.2. Velocity Dispersion Profile 

As relaxed halos are in equilibrium, the fact 
that their mass distribution develops inside-out 
automatically implies that their local velocity ten- 
sor keeps unaltered as they accrete. Indeed, cen- 
tered spheres of any arbitrary radius conserve not 
only the mass but also the kinetic energy (and the 
gravitational energy as well provided the origin of 
the potential remains unchanged; see the discus- 
sion below). Consequently, all kinematic profiles, 
in particular the velocity dispersion profile, must 
develop inside-out just like the spherically aver- 
aged density profile. Therefore, the total energy 
is given by (see eq. |A32| ) 

r R(t) 

E(t) = E[R(t)} = 4tt / drr 2 (p)[r) 
Jo 



cr 2 (r) GM(r) 



(21) 



with all the functions within the integral on the 
right independent of time and, by differentiation, 
we obtain 



<j 2 (t) = a 2 [R(t)} =2 



E GM(t) 
M + W) 



(22) 
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Caution must be paid to the fact that the two pre- 
ceding expressions presume the gravitational po- 
tential, with origin at infinity, for the system trun- 
cated at the virial radius. This simplifies notably 
the expression of the total energy as it depends on 
the inner mass distribution only, which is particu- 
larly useful when dealing with protohalos at very 
early times (see below) . Had we not adopted that 
point of view, we would have been led to the more 
general expression 

o 2 {t)^a 2 [R{t)]=2^ + {<5>)[R{t)] (23) 
M 

where ($)(?") is the spherically averaged gravita- 
tional potential (eq. |A6j ) associated to the in gen- 
eral non-truncated mass distribution. 

In any case, the same reasoning above lead- 
ing to the uniqueness of the (p) profile for halos 
with the same value of M at t accreting mass at 
the same rate M during some finite time interval 
around t leads now to the uniqueness of their a 
profile if they also accrete energy at the same rate 
E during that interval. In other words, the (an- 
alytical) a profile also permanently adapts to the 
energy accretion currently undergone by the halo. 

2.3. Anisotropy Profile 

According to the preceding discussion, the (p) 
and <7 profiles of a relaxed halo are completely 
set by the independent mass and energy accre- 
tions it is currently undergoing. What enables 
the emerging structure to be in equilibrium is 
the freedom provided by the anisotropy, which 
adapts to produce a steady configuration consis- 
tent with the inside-out growth. More specifically, 
for some given (p) and a profiles, the (approxi- 
mate) Jeans equation (|A30[) for anisotropic sys- 
tems is a differential equation for c r (r), which 
fixes the anisotropy profile. To see it consider the 
virial relation (|17[) resulting from integration of 
the Jeans equation. Taking into account the rela- 
tion (|2T1), we can write 



x 2E(t) + 



R{t) ,GM(r) 
dM(r) ^ 



(24) 



by both independent (p) and a profiles down to 
the halo center. 

To sum up, the fact that relaxed halos develop 
inside-out implies that those having identical val- 
ues of M and E at t and undergoing the respective 
accretions at the same rates M and E during any 
finite time interval around t necessarily have not 
only identical (p) and a profiles, but also identical 
er r (and, hence, at as well) and (3 profiles. 

2.4. Angular Momentum Profile 



As mentioned, the expression (|13|) for the ki- 
netic energy within r includes the influence of ro- 
tation. As the mass and kinetic energy of halos in 
centered spheres of arbitrary radii are kept unal- 
tered during accretion, the total angular momen- 
tum inside them must also necessarily be kept un- 
changed. The total angular momentum of accret- 
ing halos may, of course, change with time. But, 
according to the preceding reasoning, this can only 
be done at the expense of the matter newly incor- 
porated at the instantaneous radius like for the 
total mass and energy of the system. 

Under these circumstances, equation (TIrJ]) leads 

to 

J(t) = J[R(t)] = 4tt f m drr 2 (j)(r)(p)(r) , 
Jo 

(25) 

where all the functions inside the integral on the 
right are again independent of time and, by differ- 
entiating, we obtain the relation 



(j)(t) = {?>[#(*)] = i 



(26) 



Thus, the value of the radial velocity dispersion at 
any given radius is completely determined indeed 



Then, the same reasoning leading to the unique- 
ness of the (p), a and j3 profiles, added now to the 
approximation of aligned local angular momenta, 
leads to the fact that halos with given values of M 
and J (or of M, E and A or A') at t, increasing at 
the same respective rates M, and J (or M, E and 
A or A') during any finite time interval around t, 
also have identical (j) profile at all radii. 

3. TYPICAL PROFILES 

The fact that all the preceding profiles are 
unique for given current values of the correspond- 
ing accretion rates is not of much help, in general, 
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Fig. 1. — Predicted density profiles (top panel) 
as a function of radial distance in units of virial 
radius, R, for halos at z — (black solid lines), 
compared to their best fit by a NFW law (red 
dashed lines) and by a Sersic law (blue dotted 
lines) down to R/100 and f pc, respectively, for 
halos with 10 10 , 10 u , 10 12 , 10 13 , 10 14 and 10 15 
Mq (from bottom to top). The normalization con- 
stant jo, is defined as pM*/M where M* = 10 13 
M Q is the critical mass for collapse. Correspond- 
ing radial dependence of logarithmic slopes (bot- 
tom panel) for the same halo masses: fO 10 M 
(dot-long dashed line), 10 11 M Q (dot-short dashed 
line), 10 12 M (long dashed line), 10 13 M (short 
dashed line), fO 14 M Q (dotted line) and fO 15 M 
(solid line). 



for inferring them for any particular halo. This re- 
quires performing the analytical extension of the 
respective small pieces developed by accretion in 
any finite time interval, which is not an easy task. 
However, what can be readily derived in the way 
explained next is the typical radial behavior of all 
these profiles. 

3.1. Density Profile 

As noted by SMGH, the analytical extension of 
the typical mass profile M(r) is simply the trans- 
formation from t to r through the inverse of equa- 
tion (|2"0)) of the analytical M(t) track solution of 
the differential equation 



M 



= K a [M(t),t] 



(27) 



where 



n a [M(t),t] 



MA„ 



d(AAf) AM K m {M, t, AM) 
. (28) 

is the analytical typical scaled accretion rate in the 
cosmology under consideration (Raig et al. 2001). 
Then, one must simply differentiate such a mass 
profile to obtain the desired typical density profile, 



(p)(t) 



t {t)p(t) 
1 



dln(A vir p) 



n a [M(t),t] 



dt 



(29) 



for t{r) equal to the inverse of R(t) given by equa- 
tion ([ll- In equation (25]), TZ m (M, t, AM) is the 
usual merger rate in the extended Press-Schechter 
(PS) formalism (Lacey & Cole 1993) and the in- 
tegral extends only over those mergers producing 
a relative mass increase below the effective ratio 
A m = AM/M separating minor from major merg- 
ers as the former are the only ones to contribute 
accretion. As shown in SMGH, the best effective 
value of this threshold, equal to 0.26, gives an ex- 
cellent fit to the empirical NFW c(M) relation 
over at least 4 decades in halo mass (see Fig. 2 
in SMGH). 

In the top panel of Figure [1] we show the pre- 
dicted density profiles, for A m = 0.26, correspond- 
ing to various halo masses at z = 0, compared to 
their respective best fits down to a radius of one 
hundredth R and 1 pc by a NFW law (eq. [1]) 
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Fig. 2. — Predicted phase-space density profiles at z = for the same halo masses as in Figure [T] (and same 
symbols as in the bottom panel of that figure). The red thick long dashed-short dashed line is the empirical 
curve obtained by Ascasibar et al. (2004). Left panel: predictions obtained from the spherical collapse 
model without shell-crossing. Right panel: predictions obtained from the phenomenological correction for 
shell-crossing with r/ — 1.8. The density and velocity dispersion are in units of the cosmic critical density 
and of halo circular velocity, respectively. 



and an Einasto law (eq. [2]), respectively. As men- 
tioned, not only is the density profile well fit by an 
NFW law down to one hundredth the virial radius 
as found in numerical simulations, but the mass- 
concentration relation is also perfectly recovered. 
In respect to the fit by an Einasto law, there is 
so far no accurate mass dependence of the corre- 
sponding parameters in the literature; see SMGH 
for that predicted by the present model. In the 
bottom panel of Figure [T] we plot the correspond- 
ing logarithmic slopes as a function of radius, 
where the tendency to approach, at small radii, 
the functionality dln(p)/dlnr oc — r 1 /" character- 
istic of the Einasto law is more apparent. 

3.2. Velocity Dispersion Profile 

To derive the typical a profile we need to know 
the typical energy accretion rate E going along 
with the typical mass accretion rate M given 
above. Although there is no expression similar 
to (|2"T|) for E, we can try to estimate it by means 
of conservation arguments from the typical total 
energy of the halo seeds at any arbitrarily small 
time t\. 

When the total mass M(t) of virialized objects 
at t is to be estimated it is very useful to ap- 
proximate protohalos by smooth spherical top-hat 



perturbations. This does not work however when 
dealing with the total energy E(t) because it de- 
pends on the halo mass and velocity distributions 
at all scales. (Nor does it help taking into account 
that protohalos would coincide with peaks of the 
smoothed density field.) Determining accurately 
the total energy of the protohalo is, in general, a 
very difficult task. But, in the present case, we 
can take advantage of the fact that the structure 
and kinematics of the final halo do not depend on 
its past aggregation history and adopt the point of 
view that it has evolved since ti by pure accretion. 
That assumption and the neglect of triaxiality ef- 
fects as done in the Appendix for the final relaxed 
halofl allow one to deal with those seeds as if they 
were spherically symmetric. 

In this case, the radius R p h of that part of the 
protohalo collapsing at t is 

3 _ 3M(t) 

R ^ - 4np(m+s(tM ' (30) 

where S(t,U) — 6 c D(ti)/D(t) is the density con- 
trast for spherical collapse at t, with S c the criti- 
cal value for current collapse (equal to 1.69 in any 



2 Such accreting protohalos are also much centrally peaked, 
which guarantees the validity of the approximation. 



8 



3 

O 







II 


ii 


i i i i i l_ 












- 


III 


ii 


ii 


-2 


-1.5 -1 


-0.5 



log(r/R) 



Fig. 3. — Predicted phase-space density profile for 
a halo of 10 13 M Q at z = for different values of 77: 
(long dashed line), 1 (dot-dashed line) 1.5 (short 
dashed line), 1.8 (solid line) and 2 (dotted line) 
compared to the empirical curve obtained by As- 
casibar et al. (2004) (red thick long dashed-short 
dashed line). 



flat cosmology) and D{t) the perturbation linear 
growth factor. Consequently, its mass and total 
energy (for the system truncated at i? p h and with 
potential origin at infinity) are 



M{t) = 4tt 



■Rph(t) 



Pph(r) 



(31) 



r-Hph(t) 

E(t)=4irj drr 2 p ph (r) 



Mp h (r) GM ph {r) 



(32) 

where p P h(r) and M p h(r) are the exact, spherically 
averaged, density and mass profiles, respectively, 
of the protohalo and u p h{r) is the velocity of the 
shell at r, essentially equal to the Hubble com- 
ponent, H(ti)r. Differentiating over R p h(t) both 
M(t) and E(t), respectively given by equations 
(f3Tj) and (f32"l) , and substituting these derivatives 
into equation (|22jl . one is led, after some algebra 



and to leading order in the perturbation S(t, ij), to 
GM(t) 



a\t) = 2- 



R(t) 



1 - 



(33) 



yielding, for t equal to the inverse of R(t), the 
wanted typical <J 2 (r) profile. Equation (f33|) coin- 
cides with the expression we would have obtained 
from the usual top-hat approximation, so the de- 
tailed density and velocity distributions in the pro- 
tohalo do not actually play any role in the final 
result. The reason for this is that the specific 
density profile /0 p h(f) cancels when taking the ra- 
tio between the time derivatives of E and M in 
equation (|2"2"|) . This is the consequence of the fact 
that, as the velocity dispersion profile is develop- 
ing inside-out, its value at r depends only on the 
specific energy of the protohalo at the correspond- 
ing radius, not on its total inner value. 

The phase-space density profile (p) /a 3 arising 
from such a velocity dispersion profile and the 
spherically averaged density profile derived in Sec- 
tion 13.11 is shown in the left panel of Figure [2] 
For comparison, we plot the empirical best-fitting 
power law of index v = 1.9 obtained by Ascasibar 
et al. (2004). The theoretical prediction agrees 
with the result of numerical simulations (despite 
the lack of any free parameter to adjust) not only 
in its overall shape, close to a power law over more 
than five decades in halo mass, but also in its mag- 
nitude. The only small discrepancy is that the 
predicted profile is somewhat steeper than empir- 
ically found (v ~ 2.3 instead of 1.9). This is most 
likely due to the fact that we have neglected shell- 
crossing. 

Shell-crossing brakes the infall of collapsing 
matter, so the bulk velocity of deep enough lay- 
ers is much smaller than predicted by the top- 
hat collapse model without shell-crossing. In fact, 
the system is essentially steady within the virial 
radius R (Cole & Lacey 1996). Unfortunately, 
shell-crossing cannot be dealt with analytically, 
so the mass within R(t) can only be estimated 
phcnomenologically from the predictions of that 
simple collapse model. Next a similar approach 
is used to estimate the total energy of the object 
within R(t). 

The mass within the shell of the protohalo at 
i?ph having an inner mean density contrast ap- 
propriate for collapse, without shell-crossing, at t 
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appears to coincide with the mass M of the fi- 
nal steady object within the virial radius R, de- 
fined through equation (|20p. at that moment. This 
does not mean, of course, that the particles in- 
side R coincide with those initially located inside 
i? p h- Some have rebound and are currently be- 
yond R(t), while others initially beyond i? p h have 
already passed by R(t). Yet, as far as the mass 
within R(t) coincides with that of the protohalo 
within i? p h, the mass of these two kinds of particles 
balance each other and we must not worry about 
that distinction. However, the situation is differ- 
ent when dealing with total energy. Particles hav- 
ing bounced were originally more tightly bound 
than those having not yet. Thus, the total en- 
ergy within R(t) is larger than estimated through 
equation (|32p (kinetic energy transfer among shells 
through two-body interactions at shell crossing is 
negligible). In other words, it is given by that 
equation at a larger time and the same result holds 
for the velocity dispersion at R: it is given by equa- 
tion (|33|) at some later time t + At. Inspired of the 
fact that the shift At in the case of M(t) is null 
for all halo masses and cosmologies, a reasonable 
guess for the shift At in the case of energies is to 
assume it equal to rjt with r\ equal to the same 
positive constant factor for all halo masses and 
cosmologies. This leads to (see eq. [33] ) 



a 2 (t) = 2 



GM[(l + ri)t] 
R[(l + v )t] 



1/3 



[(!+»»)*] 



(34) 

Moreover, since the shift rjt is caused by shell- 
crossing and the characteristic bounce time of 
shells that would collapse at t if there were no 
shell-crossing is equal to t, we expect an rj value 
of order unity. 

In the right panel of Figure [2J we depict the 
phase-space density profile arising from the spa- 
tial velocity dispersion given in equation (|34[) for 
t(r) equal to the inverse of R(t) and rj = 1.8. As 
can be seen, this new a profile is in much better 
agreement, indeed, with the results of simulations. 

To see the effect of varying 77, we show, in Fig- 
ure [31 the profiles resulting from different values 
of that parameter. For any value of 77 greater then 
zero, the solution <r(r) at large enough radii re- 
lies on the extrapolation towards the future of the 
Bryan and Norman (1998) expression for A v ; r . For 
rj = 1.8, this affects the solutions at radii above 



log(r/i?) ~ —0.35 for all halo masses (see Figs. [2] 
and[3]where all profiles are plotted only for smaller 
radii). As such an extrapolation is quite uncertain, 
from now on, the corrected a profile at such large 
radii is inferred from the much more secure linear 
extrapolation of the predicted phase-space density 
profile. 

3.3. Anisotropy Profile 

After a change of integration variable, the ap- 
proximate relation |24[) takes the form 



x / dt 
la 



1 



A^{t) {p){t) 



GM(i) 



R{t) 



n a [M(t),t\M(t), (35) 



leading, for t(r) equal to the inverse of R(t), to 
the (approximate) typical radial velocity disper- 
sion profile and, through equation [[%]). to the (ap- 
proximate) typical P(r) profile. Note that, as a 
was derived to leading order in the perturbation 5 
at tj, its central behavior is not fully reliable, so is 
not cither that of the c r profile. However, the be- 
havior at moderate and large radii of <7 r (r) is quite 
robust because the integral appearing in equation 
([35| is very insensitive to the central values of a. 

In the left panel of Figure [4] we plot the pre- 
dicted typical (3(r) profile for the same current ha- 
los and the same radii as in Figure [5] In the right 
panel, we show the corresponding relation between 
(i and the logarithmic slope of the density profile. 
The theoretical anisotropy shows a trend similar to 
that empirically found by Hansen & Moore (2006), 
although with some apparent undulations. Similar 
variations have recently been observed in numeri- 
cal experiments (McMillan et al. 2007). Whether 
these undulations are real or an artifact due to 
the approximate character of our solution is hard 
to tell. Note that the rough universality of the j3 
profile combined with the fact that the typical halo 
density profile at t depends only on the halo mass 
leads to the conclusion that the typical anisotropy 
profile depends essentially only on the halo mass, 
too. 

Of course, such a predicted typical anisotropy 
profile depends on the exact value of 77 chosen to 
correct the theoretical velocity dispersion profile 
for the effects of shell-crossing. To see more quan- 
titatively the influence of such a correction, we 
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Fig. 4. — Predicted anisotropy profile as a function of radius (left panel) or density logarithmic slope (right 
panel) obtained for the same halo masses as in Figure [2] For comparison, we plot the best linear fit to 
empirical data on this latter relation drawn from numerical simulations (central red thick short dashed-long 
dashed line) and its scatter (bracketing red thick lines). For r tending to zero, the predicted anisotropy tends 
to become negative. Whether this is a real effect or a numerical artifact (caused, e.g., by the extrapolation 
of our phenomenological correction for shell crossing down to r — 0) is hard to tell at this stage. 




Fig. 5. — Same as Figure [3] for the predicted anisotropy profile, showing the effect of the phenomenological 
factor r/. Full line (77 = 1.8) gave a good fit when comparing the theoretical (p)/cr 3 to those of simulations, 
and it is also seen to provide a reasonable fit to (3. 
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Fig. 6. — Top panel: predicted phase-space densi- 
ties for a halo of 10 13 M© at z — correspond- 
ing to the radial (er x = <r r , dashed line), tan- 
gential (<r x = at, dot-dashed line) and spatial 
(er x = a, solid line) velocity dispersions obtained 
with r\ — 1.8 as in the right panel of Figure [U 
Bottom panel: residuals from the respective fit to 
a power-law of the form Ar~ v with best-fitting ex- 
ponents equal to 1.89 (spatial case), 1.92 (radial 
case) and 1.87 (tangential case); all fits achieved 
in the radial range requiring no extrapolation. 



Fig. 7. — Mass dependence of the power index v 
(top panel) and the proportionality factor A (bot- 
tom panel) of the best-fitting power law Ar~ v for 
the predicted spatial, radial and tangential phase- 
space densities (same symbols as in Fig. [6]). 
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plot in Figure [5] the solutions arising from the dif- 
ferent values of rj used in Figure [3] 

Once we have determined the typical radial and 
tangential velocity dispersion profiles we can check 
if the corresponding phase-space density profiles 
are again close to power laws. The result is shown 
in Figure [6] where we also plot the residuals of each 
profile from the corresponding best fit by a power- 
law in r (eq. [3] ) . As can be seen, the phase-space 
density profiles for the radial and tangential ve- 
locity dispersions also admit a power law fit like 
that associated with the spatial velocity dispersion 
with similar values of index v and just a small 
shift in the respective proportionality factors, in 
agreement with the results of numerical simula- 
tions (Ascasibar et al 2004; Dehnen & McLaugh- 
lin 2005). Among these two profiles that associ- 
ated with the tangential dispersion, (p)/af , is the 
one giving the best fit to a power-law in radius 
(in the range requiring no extrapolation), with a 
slope v ~ 1.9. This result is most likely related to 
the fact that the tangential velocity distribution 
function has virtually the same shape for all radii 
(Hansen et al. 2006). This is contrasted with the 
radial velocity distribution function, whose shape 
changes significantly as function of radius. 

In Figure [7] we show the mass dependence of 
the best fitting parameters v and A to all three 
phase-space density profiles. 

3.4. Angular Momentum Profile 

Equations (f26|) and (|27|) define, for t(r) equal 
to the inverse of R(t), the (j)(r) profile in terms 
of J, which can be computed from equation ^ 
or ||7J), assuming typical values for both the spin 
parameter and its time derivative. 

As mentioned, the mean spin values and stan- 
dard deviations are found in iV-body simulations 
to be quite insensitive to the cosmology, halo mass 
and particular epoch considered. When looking at 
their behavior in more detail, it is observed how- 
ever that, just after major mergers, they take val- 
ues substantially higher than the mean (Burkert 
& D'Onghia 2004), which is likely due to the fact 
that halos are not yet fully relaxed. What is more 
important for our purposes here, A' depends on 
halo mass according to equation |[H} (BOlb), its 
mean value being essentially constant while the 
mean A value is essentially independent of mass 



(BOlb) and shows a slight secular evolution (Het- 
znecker and Burkert 2006). This latter result is 
found regardless of whether the means are per- 
formed over all the halos or just accreting ones. 
According to the relation (jHJ), the spin parameters 
of a given accreting halo satisfy the relation 



X'(t) = X(t)F{c[M(t),t]}. 



(36) 



Taking the mean of this latter equation over halo 
masses in the relevant mass range at t and tak- 
ing into account that the mean A' is constant and 
equal to 0.039 one is led to 



A(t) = 



0.039 

~W) 



(37) 



where F(t) stands for the mean of F[c(M, t)] with 
F(c) given by equation © and c(M, t) the NFW 
mass-concentration relation at t given in SMGH 
(see their eq. [10]). Note that, at present, c is 
around 10 for masses in the relevant range (be- 
tween 10 10 M Q and 10 15 M ), so F is about unity. 
We have checked that the time-dependence of the 
mean A value implied by the approximate relations 
(f3"6"]) - (|3"T)) is also a reasonable approximation of the 
one found by Hetznecker and Burkert (2006). 

Differentiating equations ([6]) and 0, we obtain 
the typical angular momentum accretion rate 



J 

W) 

and 
j 

W) 



M \b 



1 



M{t) 1 2 



M \b 

M(t) j 3 



E 



n a [M(t),t] 



1 

dln(A vir p) 1 / 6 
df 



2E{t) 



A 

W) 



(38) 



X'(t) 



(39) 



respectively. In equation l|38[). the energy accre- 
tion rate E is computed from equation ([23]), with 
a(t) according to the expression (f3"4"| and ($) the 
spherically averaged gravitational potential at R 
of the halo endowed with a NFW profile up to in- 
finity, and E(t) is obtained by integrating E over 



3 Strictly, this value corresponds to the mean for all halos, 
that for accreting halos being slightly smaller (Hetznecker 
and Burkert 2006). However, this is a very small effect that 
can be neglected at this stage. 
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Fig. 8. — Predicted specific angular momentum profile for the same halo masses (marked with identical 
symbols) as in Figures [2] and |4] derived using the A (left panel) and A' (right panel) spin parameters. For 
comparison, we plot the average power law with s = 1.3 fitting, to a first approximation, numerical data 
according to BOlb (red thick short dashed-long dashed line). 



time. In deriving equation (|39p from equation ((T)), 
we have taken into account the definition of the 
virial radius (eq. [20] ) and the expression of (p) (t) 
given by equation (f29|) . Note that the expression 
(f3"5)) making use of A' does not depend neither on 
the phcnomcnological correction for shell-crossing 
of the a profile nor on the assumed mass distribu- 
tion beyond R(t). 

By substituting into equation (|26|) the previ- 
ous expressions for J and taking into account the 
typical values of the spin parameters given by ex- 
pressions ([36]) and ([37]) . we are finally led to the 
two following, in principle equivalent, estimates of 
the typical spherically symmetric specific angular 
momentum profile, 



(j)(t) 



and 



<i>(*) 



M(t)[2 K a [M(t),t] 



At) (5 
M{t) \ 3 



E 



2E(t) 



A 

W) 



(40) 



n a [M{t),t] 

dln(A vir p) 1 / 6 
dt 



A'(t) 



(41) 



for t(r) equal to the inverse of R(t) and J(t) the 
total angular momentum accretion track, given by 
equations ^ or ([7]) for the appropriate values of 
M(t), E(t) and A(t) or X'(t). 

These two estimates of the angular momentum 
profile are shown in Figure [8] As can be seen, 
the solutions derived from A and A' are slightly 
different from each other. This cannot be due to 
the value of -q used to infer the solution from A 
because, as shown in Figure [TtJI the solutions ap- 
pear to be remarkably insensitive, in this occa- 
sion, to the value of that parameter. We have also 
checked that the extrapolation of the density pro- 
file beyond R has a very small effect. Thus, it can 
only be due to the time dependence adopted for 
A and A'. In fact, the information we have on the 
empirical behavior of the spin parameters refers 
only to redshifts smaller than 2 (see Hetznecker 
and Burkert 2006). For this reason, in Figure [SJ 
the predicted solutions are only drawn for radii 
involving z < 5. 

Anyhow, both solutions show similar trends in 
agreement with the results of numerical simula- 
tions. They are indeed power laws in mass inside r 
to a first approximation, particularly the solution 
obtained from A, with a mean best-fitting value 
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Fig. 9. — Same predicted specific angular momentum profiles (and same symbols) as in Figure [3 and the 
corresponding fits either by a power law for the A case (left panel), or by the slight modification of it given 
in the text for the A' case (right panel). To avoid overcrowding, we only represent the curves corresponding 
to 10 10 M B and 10 15 M m . 



of index s (with the restriction (j)(R)/j max = 1) 
close to that found empirically. More specifically, 
for the masses included in our figures and a PS 
mass distribution, we obtain (within the radial 
range plotted in Fig. [8] for each mass) s ~ 1.6 and 
~ 1.7 for the A and A' cases, respectively. While 
for the mass interval studied by BOlb (> 10 12 
Mq) we obtain s ~ 1.3 and ~ 1.4 to be com- 
pared with the value of ~ 1.3 quoted by these au- 
thors. Furthermore, the slight bending upwards 
at large radii, more apparent in the A' case al- 
though also present in the A one as all the curves 
must approach unity at r = R, is also observed in 
simulated halos (see the comments in BOlb). Fi- 
nally, the tendency for less massive halos to have 
slightly steeper profiles (in the mean all over the 
radial range analyzed) is also consistent with the 
empirical results of BOlb (see their Fig. 16). 

As a power law does not seem to give a very 
good fit to the solutions drawn from A' we have 
tried with another simple analytical expression. 
As shown in Figure [H the following slight mod- 
ification 



(Mr) 



= 1 - 



1 - 



M{r) 
M(R) 



(42) 



is already enough to improve notably the fits for all 
halo masses. Note that the solutions drawn from A 
are however better fit by the previous simple power 
law of index s. The dependence on M of that index 
s and the new one s' is given in Figure [TT] for the 
solutions drawn from A and A', respectively. 

4. SUMMARY AND CONCLUSIONS 

We consider the structural and kinematic prop- 
erties of realistic ellipsoidal, rotating, standard 
CDM halos in the accretion-driven scenario. Ac- 
cording to this scenario, halos adapt dynamically 
to the amount of matter (carrying mass, energy 
and angular momentum) that is accreted at any 
given moment. This is due to the collisionless na- 
ture of CDM and the inside-out growth of these 
structures during accretion phases. The typical 
accretion rate for halos of a given mass at a given 
cosmic time depends only on the cosmology, and 
thus, the typical radial profiles of relaxed dark 
matter halos are fixed by cosmological parameters. 

Specifically, we derive the typical spherically 
averaged profiles corresponding to the density p, 
the velocity dispersion a and the phase-space den- 
sity p/cr 3 . We find that all these predicted profiles 
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Fig. 10. — Same as Figures O and [3] but for the 
angular momentum profile predicted using A for 
the extreme cases of 10 10 Mq (lower curves) and 
10 15 M Q (upper curves) in order to show the de- 
gree of insensitivity to parameter 77 for solutions 
corresponding to different halo masses. 
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Fig. 11. — Mass dependence of the best fitting pa- 
rameters s and s' in the power law and the mod- 
ification of it proposed to fit the specific angular 
momentum profile drawn from A (solid line) and 
A' (dashed line), respectively. 



are in good agreement with those found in numer- 
ical simulations. In particular, p is well fit, in 3D, 
by an Einasto function and p/<7 3 by a power law, 
Ar~ v , with v rj 1.9 and \o%A ps 1.4. 

We also derive the typical profile of the veloc- 
ity anisotropy (3. We find that the anisotropy in- 
creases from something small (maybe slightly neg- 
ative) at the halo center to a positive value in 
the outer region and with a typical logarithmic 
slope also in agreement with numerical data. This 
is the first time that a non-zero anisotropy has 
been derived from first principles. We find that 
the spherically averaged phase-space density asso- 
ciated to the tangential (radial) velocity disper- 
sion, p/af (p/af) is also close to a power law with 
slope v w 1.9 although somewhat different values 
of A. 

We finally derive the typical spherically aver- 
aged profile of the specific angular momentum j 
from the empirical evolution of the dimensionless 
halo spin parameter A or A'. We find that this pro- 
file scales, to a first approximation, as a power law 
in inner mass, j(r) oc [Af(r)] s , in agreement with 
numerical simulations. The average logarithmic 
slope s found using A coincides with the empirical 
value of 1.3, while that found using A' (1.4) is also 
close to it. A more accurate analytical expression 
is provided for the j(r) profile derived from A'. 

We find that all the preceding profiles depend 
moderately on halo mass. The dependence of the 
NFW and Einasto parameters fitting the density 
profile in any given cosmology were already pre- 
sented in SMGH. Here we have focused on the 
mass dependence of the remaining profiles. The 
behavior of any of these profiles found in ob- 
servations (for instance, the logarithmic slope of 
the density profile obtained from X-rays or strong 
gravitational lensing) could in principle be used, 
in a given cosmology, to estimate the total mass 
of the system. Since the total mass most often is 
known, this also provides a direct way of obser- 
vationally testing the predictions of this accretion 
driven scenario. 

The anisotropy and velocity dispersion (or 
phase-space density) profiles are derived using a 
practical phcnomenological correction for the ef- 
fects of the complex process of shell-crossing on 
the total energy of halos conserved since the epoch 
they were density perturbations. This correction 
is achieved by adjusting the value of one free pa- 
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rameter, r\. We find that both the anisotropy 
profile and the phase-space density profile are in 
good agreement with numerical results for r/ w 1.8, 
the angular momentum profile being instead in- 
sensitive to the value of this parameter. 

The fact that all our results are in fairly good 
agreement with those drawn from numerical cos- 
mological simulations leads support to the idea 
that the structural and kinematic properties of 
dark matter structures are determined by their dy- 
namical adaption to cosmological accretion. 
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A. JEANS EQUATION AND SCALAR VIRIAL RELATION 
A.l. Exact Relations 

The internal dynamics of relaxed inside-out evolving halos should be well-described by a steady (true) 
phase-space density /(r, v) satisfying the collisionless Boltzmann equation. Writing this equation in spherical 
coordinates (see, e.g., equation [4p-2] of Binney & Tremaine 1987), multiplying it by the radial velocity v T , 
and integrating over velocity and solid angle, one is led to the following first order differential equation 

d« P K 2 ) , (p)(r) /Q _ 2 ^ 1 



dr 



r 



(3<j 2 (r) - cr 2 (r)) = -— / dtp &6 sin6> pd r $ , (Al) 
47r Jo Jo 



where a stands for the non-centered velocity dispersion in the shell at r, 



1 C27T 



4tt(p) 



ii 



a 2 (r) = — -- / d^l d9 sm9 I dW / = <7 2 (r) + of (r) + a 2 (r) , (A2) 



(p) for the spherically averaged density profile of the halo given by equation (|12p. with the local density 
distribution satisfying the relation 

p(r,9,<p) = J d 3 V /(r,v), (A3) 



and d t for the radial partial derivative. To derive equation (|Aip . it is only needed such conventional assump- 
tions as the continuity, in real space, of the local density and mean velocities and the fact that / vanishes for 
large velocities. In expression (|Alj) . we have taken into account that the effects of a non-null cosmological 
constant are negligible on halo scales. 

Splitting the local density and gravitational potential as 

p(r,9,<p) = (p)(r) + 5p(r,9,<p) (A4) 
<Z>(r,e,<p) = ($)(r)+6$(r,d,<p), (A5) 

and taking into account that the spherically averaged gravitational potential, 

($)(r) = — / dtp / d0 sin 6» $(r,9,<p), (A6) 
4tt Jo Jo 

satisfies, by the Gauss theorem, the usual relation for spherically symmetric systems 

d($)(r) _ GM(r) 
dr r 2 ' 

equation (|Aip adopts the form 

^1 + ,2) = - {p) ^l _ J_ f\ f\ e sin9 Sp di m . (A8) 

dr r r z 4tt J J 

Except for the second term on the right, equation (|A8|I looks exactly as the classical Jeans equation 
for spherically symmetric systems with anisotropic velocity tensor. Thus, multiplying (|A8|) by 47rr 3 dr and 
integrating over r, the same steps leading to the scalar virial relation for a spherically symmetric system now 
lead to 

AttR 3 P(R) - 2K = - f dM(r) GM ^ - f drr 3 /" dcp [ d9 sin 9 6pd T (5$) , (A9) 
Jo r Jo Jo Jo 

where dM(r) = 47rr 2 (p)(r)dr is the mass of the elementary spherical shell of radius r, 

P(R) = (p)(R)a?(R) (A10) 



(A7) 
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is the spherically averaged radial boundary pressure, and 



K = l 
2 



is the kinetic energy within R. 

On the other hand, the total potential energy is 

r-R 



dM(r)a 2 (r) 



(All) 



W=- drr 2 dp / d6 sm9 p(r,6,<p) <f>(r,6,<p) 
2 Jo Jo Jo 



(A12) 



which, from equations (IA4l) - (IA5j) . can be written as 



W = 2ir / drr 2 (p){r) ($)(r) + 



dr r 2 



dip I d6 sin0 5p8<& . 



(A13) 



Provided the central asymptotic logarithmic slope of (p) is greater than —5/2, fixing the origin of the potential 
$ so to have the boundary condition 

mR) = ~, (Au) 

for the differential equation (|A7[) (this is essentially equivalent to consider the potential origin at infinity for 
the system truncated at R) and integrating by parts (two consecutive times) the first term on the right of 
equation (|A13|1 . the potential energy takes the form 



W 



R % GM ( r ) 
dM (r) — 



dr r 2 



2tt 



dip / d6 sin 6 Sp 5$ . 



(A15) 



Therefore, by subtracting 2W on both sides of equation (|A9j) we arrive to the virial relation 



4ttR 3 P(R)-2E= [ dM(r) GM ^ _ f drr 2 /" dip f d6 sinf Sp [5$ + rd t (6$)] . (A16) 
Jo r Jo Jo Jo 



where the total energy E = K + W is given by 



E = 



dM(r) 



a 2 GM (r) 



dr r 2 



2tt 



dp / d6 sxa.6 5p5$ . 



(A17) 



A. 2. Approximate Relations 



Although halos exhibit substantial triaxiality (the average minor to major axial ratio takes a value between 
0.6 and 0.7; Bullock 2002, Kasun & Evrard 2005, Bailin & Steinmetz 2005, Libeskind et al. 2005), their mass 
distribution is far from flattened and the isopotential surfaces are more spherical (cf. Binney & Tremaine 
1987). We therefore have 



N 
(P) 

m 



< i 

< i. 



(A18) 
(A19) 



Moreover, since the isopotential surfaces approach to spheres as r increases, the rms value of (5$ in the 
spherical shell at r, 

r 2-K PTT n V 2 



— J dp J dO sin6» (8>$>f 



(A20) 
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is a positive monotonously decreasing function of r satisfying the inequality 

0"<I> 



Then, the fact that |($)| is also a monotonous positive decreasing function of r (see eqs. [A14] and |A7| ) 



dcr$(r) 



dr 



(A22) 



<1. (A21) 

he fact that \{&)\ is also a monotonous positi 

implies 

To see it, consider the proportionality 

a<s>{R) = e| ($)(#) | (A23) 

with e much smaller than unity, as implied by the inequality (|A21j) . If we move to the halo center, the positive 
increments Atr$ = u$(0) — er$(i?) and A|($)| = |($)(0)| — cannot satisfy the e-order relation 

A<7*~e"A|($)|, (A24) 

with [i any negative integer or null number, because one would then have 

cr$(0) _ cr$(i?) + A<7$ 6 + 6 + eM <J)\i] 



|(*)(0)| |($)(i?)| + A|($)| 1 I A IWI 1 I A IWI 



> 1 > (A25) 



with the last inequality arising from the condition (| A24|) and the fact that A|($)|/|($)(i?)| is of order unity 
(see eqs. |A7| and |A14| ~), which contradicts the condition (|A21[) . Thus, they must necessarily satisfy instead 

Aa$ ~ e p A|($)| (A26) 

with \i some positive integer number, in which case the relations 

da $ Aa„ A|($)| d|<*>| , A97 , 

validate the inequality (|A22|) . 

For fixed values of 9 and cj>, \6<&\ is a continuous function of r with typical amplitude equal to er$ except 
in those directions where it vanishes (this necessarily happens in specific radial directions owing to the 
homologous character of triaxial halos). We therefore have 

dan 



\d t s$\ = \d r m\ 

so equations (|A18|) and (|A22j) lead to 



dr 



(A28) 



|^||9 r 5$|«(p)(r)-^i ) (A29) 

the latter inequality also holding now over those radial directions where vanishes and, hence, d^S® is 
identically null. 

The inequality (|A29[) allows one to write equation (|A8|) in the approximate form 

d((lK 2 ) , (p)(r) ro 2 . , 2 , GM(r) 

— [3<7 r (r) - cr (r)] w -<p)(r) — 3 — (A30) 



dr 



This expression coincides with the Jeans equation for spherically symmetric systems in terms of the non- 
centered velocity dispersion profiles. Therefore, proceeding in the same way as leading to equation (|A16|) we 
obtain approximately the usual scalar virial relation for spherically symmetric systems 

4nR 3 P(R)-2E^ ( R dM{r) GM ^ , (A31) 
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where P(R) is given by equation (|A10p and the total energy E takes, from equations (|A1 1|) and (|A15p and 
the inequalities (|A18[1 - (| A19)) . approximately the same form as in the spherically symmetric case, 



E = K + W « 4tt 



drr 2 (p)(r) 



a 2 (r) GM{r) 



(A32) 
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